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Abstract — this paper introduces the application of wave variable 
method in reducing the effect of delay during the transmission of 
control signals via large distances in a multimachine power 
system. Fast response and pure control signals applied to a power 
system have great importance in control process. At first we 
design a centralized optimal control considering no transport 
delay and then we compare them to the one considering this effect 
in order to demonstrate the deteriorating effect of transmission 
delay. Then we apply the wave variable method and it is shown 
that performance of centralized optimal controller is improved. 

Keywords- Centralized optimal control, delay; multimachine; 
wave variable 



I. 



Introduction 



Over the past years, there have been many researches about 
dynamic stability of power system and many solutions have 
been introduced. The Linear Quadratic Regulator (LQR) 
optimal feedback is one of many tools to improve the stability 
of an interconnected power system. Using LQR theory, it has 
been established that for a controllable LTI system, a set of 
Power System Plant optimal feedback gains may be found 
which minimizes a quadratic index and makes the closed-loop 
system stable [1,2,3]. 

Applying the control signals to a multimachine power 
system requires to sending these signals through large distances 
of transmission channels which impose time delay to the control 
signals. Control signals should be applied fast and accurate in 
order to have their ameliorative effects. This issue is very 
important particularly during the faults and disturbances. 

Wave variable method applied in this paper significantly 
improve the performance of the optimal controller during the 
changes occur because of the transmission delay. Therefore, this 
idea gives a new practical way in stabilization large scale power 
systems via a centralized optimal control. 



II. 



Power system plant 



A power system model can be regarded as an interconnected 
system that consists of generations, transmissions and loads [1]. 



In this paper we proposed a 4 buses 2 generators system which 
can be a common multimachine plant [10]. 

The single line diagram of the multimachine model is 
described in Figure 1 . 
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Fig. 1 : Two-machine four-bus system 

III. LQR Optimal Feedback 

Here we explain briefly the LQR method used in this paper 
[4,6]. In general we should have the state space model of the 
system. A power system model can be written in state space 
equation as follows: 

X(t)=Ax(t) + Bu(t) (1) 

Y(t) = Cx(t) (2) 

Where A, B and C are system matrices and 



x = [x,,x 2 ,...,x n ] T ,u=[u 1 ,u 2 ,...,u n ] T 



The vectors x and u are considered as measurable state 
variables and input control signals. For example, assuming an 
n-machine model, that is written in variables x and u as 
follows: 

x = [AE^AE^AS^Ac^AT^ACV^ASR^AE^AV^AVs;] 1 
u=[U ei ,U gi ] T ,i = l,2,...,n 

A. Gain Determination 

Equation (1) and Equation (2) can be developed to solve an 
LQR optimal feedback solution. The Index value of LQR to be 
minimized can be written as follows: 
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J = ^j[x T (t)Qx(t) + u T (t)Ru(t)]dt (3) 



The solution of Equation (3) can be given as follows: 



S = A T S + SA-SBR _1 B T S + Q (4) 



If a time -varying positive symmetric matrix S(t) converges 
at t f = oo the solution of Equation (4) can be obtained from an 
algebraic matrix Riccati equation as follows: 



= A T S + SA-SBR -1 B T S + Q 



The gain K can be written as follows: 



K = R _1 B T S 
u = -Kx 



(5) 



(6) 



Then Equation (1) can be written as a closed loop system so 
that we can discuss the converging behavior using the 
following equation: 



x = (A-BK)x 



(7) 



B. Weighting Matrices Q and R 

Weighting matrices Q and R are important components of an 
LQR optimization process. The compositions of Q and R 
elements have great influences on system performance. Q and 
R are assumed to be a semi-positive definite matrix, and a 
positive definite matrix, respectively [4]. 

IV. Control Structure and System Model 
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Fig. 2: Interconnected system using LQR 

The power system assumed is constructed with 
interconnected two machines controlled by a LQR optimal 
controller as shown in Figure 2. The multimachine system is 
represented in classical model based on Kundur [1]. The inputs 
to the Generators are provided from the controller outputs, and 
the inputs to the centralized controller are variables of the 
generator plants. Each generator has two input signals 



u = [U e ,U,] T 



and 



ten 



output 



signals 



x = [AE^,AE^,A6,Aco,AT m ,ACV,ASR,AE f ,AV R ,AV s ] T .Then 

the interconnected system has 4 input signals and 20 output 
signals. The input nodes of the control signals are depicted in 
Figure 3 and Figure 4 for the speed-governing system and the 
excitation system. Both figures also describe the models used 
in this paper [1,10,11]. 

The turbine model [1] also considered as a simple single 
degree model suitable for this study. The whole system 
structure is illustrated in Figure 5. All parameters are given in 
TABLE I. 
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Fig. 3: Model of governing system 
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Fig. 4: Model of excitation system 

V. Centralized Optimal Control (COC) Results 

A. COC Applied With No Delay 

The results are in response of a symmetrical three phase 
fault applied to the Bus 4 of the proposed system during 0.5 s 
considering breaker operation time and system restoration. 
System base is given in 1000 MVA. 

In this part we assume that the delay in transmitting the 
control signals is negligible so the delay is not considered in 
both sending and receiving signals. 

From the Figures 6-9, we can see the power angle and 
voltage variations. The Optimal control response is obviously 
the best achieved. Both voltage and frequency overshoots 
damp to the steady state sooner so the settling times improve. 
The control signals applied to the machines are shown in 
Figures 10-11. 

The optimal controller has an obvious advantage in 
stabilizing the power system. 
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Fig. 5: The complete system structure 



B. COC Applied Considering Transmission Delay 

In this part the same conditions above are assumed. The 
LQR optimal controller is a centre which is positioned in some 
distance of the synchronous generators connected to centralized 
controller in its area. Sending and receiving signals (as 
mentioned before) through these large distances in a large scale 
power system imposes the transport delay which can make a 
change to the control signals and also the measured state 
variables achieved from the generators. This change 
deteriorates the stabilizing effect of the optimal control and 
seriously damages the stability of the whole system. The 
control signals are applied to the two major control systems 
which are the governing system and the AVR system as 
described before. 

From the Figures 12-15, we can see the deteriorating effect 
of transmission delay both in power angle and voltage stability 
cases. The optimal controller response is very bad and both 
generators are driven into unstable conditions. 

The voltage responses in both generators are oscillatory 
unstable and the power angles of each generators have 
increasing manner which means that the power angle stability 
is lost. The control signals applied to the machines are shown 
in Figure. 16. 

The delay time considered here is 2.22 up to 2.25 ms 
approximately. The delay time is considered to be constant 
during the simulation. We assumed that the distance between 
the machines and the centralized optimal controller is the same 
which is only for convenience and has no effect on the whole 
concept. 



VI. Wave variables method 



One of the useful methods for delay reduction is wave 
variables method. In this method, by locating two transform 
functions on two sides of the signal path (between master and 
slave), delay is reduced. Wave variables present a modification 
or extension to the theory of passivity which creates robustness 
to arbitrary time delays. The theory and its calculations are 
fully described in [5,9]. 

Here we introduce the wave variable transform. The 
parameter u is defined as the forward or right moving wave, 
while v is defined as the backward or left moving wave. The 
wave variables (u,v) can be computed from the standard power 
variables (x, F) by the following transformation. 



bx + F 

II - , ,V: 



bx-F 



(8) 



Where b is the characteristic wave impedance and may be a 
positive constant or a symmetric positive definite matrix. It 
also has the role of a tuning parameter. We also have: 



F = bx-V2bv,u = -v + V2bx 



(9) 



All other combinations are also possible. Notice especially 
the apparent damping element of magnitude b within the 
transformation. It retains the power necessary to compose the 
wave signal and then implicitly transmits the energy with the 
wave to the remote location [5]. 
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Fig. 6: Power angle variations of Gl 




Fig. 7: Voltage variations of Gl 
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Fig. 8: Power angle variations of G2 
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Fig. 12: The effect of delay on the power angle of Gl 
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Fig. 13: The effect of delay on the terminal voltage of Gl 
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Fig. 14: The effect of delay on the power angle of G2 
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Fig. 15: The effect of delay on the terminal voltage of G2 



Fig. 9: Voltage variations of G2 
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Fig. 10: Optimal control signal applied to Gl 




Fig. 1 1 : Optimal control signal applied to G2 
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control signals that should be applied to generators. This 
process happens simultaneously as the control centre always 
monitors all the variations of the generators. 

Considering wave variables transform, the right moving and 
the left moving signals should have the same dimensions. For 
this reason we designed a kalman state estimator [12] in the 
optimal control centre that by using the generators outputs as 
measurement signals and some predetermined inputs of the 
generators as known inputs, estimates the state variables 
necessary for the LQR control. Notice that the estimator is 
designed fast enough to estimate the state variables correctly. 
Now we have two control inputs as right moving signals and 
two outputs from each generator as left moving signals. 

The same fault condition assumed before is applied to the 
system and the delay time is considered 0.1 s and the parameter 
b is equal to 50. The Figures 18-21 show the results. 

2.5 r 




Fig. 18: Power angle variations using wave variables in Gl 




Fig. 16: Optimal control signal applied to Gland G2 
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Fig. 19: Voltage variations using wave variables in Gl 
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VII. Control strategy using wave variables 

Figure 17 depicts the wave variable method used in this 
paper. The master side which is the optimal control centre 
receives the outputs of the generator which are the active 
power P and the reactive power Q variations of the each 
generator. Then the optimal control centre determines the 
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Fig 17: Wave variable structure 
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Fig. 20: Power angle variations using wave variables in G2 
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Fig. 2 1 : Voltage variations using wave variables in G2 



TABLE I: Machines Data 



20 



Parameter 


Generator 1 


Generator 2 


Type 


Steam 


Steam 


Rated MVA 


920.35 


911 


Rated KV 


18 


26 


Power Factor 


0.9 


0.9 


Speed (rpm) 


3600 


3600 


D(pu) 


2 


2 


H(s) 


3.77 


2.5 


r s (P u ) 


0.0048 


0.001 


X b (pu) 


0.215 


0.154 


x d (pu) 


1.79 


2.04 


x' d (pu) 


0.355 


0.266 


Xq (PU) 


1.66 


1.96 


x'c (PU) 


0.570 


0.262 


K A (pu) 


50 


50 


T A (s) 


0.07 


0.06 


vr x (pu) 


1 


1 


V R min (Pu) 


-1 


-1 


K E (pu) 


-0.0465 


-0.0393 


T E (s) 


0.052 


0.44 


A Ex (pu) 


0.0012 


0.0013 


B Ex (Pu) 


1.264 


1.156 


K f (pu) 


0.0832 


0.07 


T, (s) 


1 


1 


Tdo (s) 


7.9 


6 


Too (s) 


0.41 


0.9 


K G 


20 


20 


T SR (s) 


0.1 


0.1 


T SM (S) 


0.3 


0.3 


T tu (s) 


0.3 


0.27 


TABLE II: Impedances of tie lines and admittance of load 


Z 14 (pu) 


0.004 +J0.1 


Z 24 (pu) 


0.004 +J0.1 


Z 34 (pu) 


0.008 +J0.3 


y 40 (pu) (load bus) 


1.2-J0.6 



VIII. CONCLUSION 

In this paper we discussed a new method in power system 
stabilization. We have considered the real situation of large 
scale power systems considering the deteriorating effect of 
time delay in the method proposed. Wave variables method 
proposed in this paper has the obvious ameliorating effect and 
reduces the effect of delay. 
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